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Abstract: To leading order in perturbation theory, we solve QCD, defined on a 
small three sphere in the large N and Nf limit, at finite chemical potential and 
map out the phase diagram in the (//, T) plane. The action of QCD is complex 
in the presence of a non-zero quark chemical potential which results in the sign 
problem for lattice simulations. In the large N theory, which at low temperatures 
becomes a conventional unitary matrix model with a complex action, we find that 
the dominant contribution to the functional integral comes from complexified gauge 
field configurations. For this reason the eigenvalues of the Polyakov line lie off the 
unit circle on a contour in the complex plane. We find at low temperatures that as fi 
passes one of the quark energy levels there is a third-order Gross- Witten transition 
from a confined to a deconfined phase and back again giving rise to a rich phase 
structure. We compare a range of physical observables in the large N theory to 
those calculated numerically in the theory with N = 3. In the latter case there are 
no genuine phase transitions in a finite volume but nevertheless the observables are 
remarkably similar to the large N theory. 
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1. Introduction 



The phase diagram of QCD at finite quark number density is of considerable interest 
but a first principles derivation of the grand potential has eluded us for two very 
good reasons. The first is that the phase transitions are conjectured to occur at 
densities where the coupling strength of QCD is large and thus the application of 
conventional perturbation theory is not valid. The other reason is that the action of 
QCD is complex resulting in the sign problem. This prevents the use of the usual 
technique of importance sampling in lattice gauge theory simulations as it is not 
possible to formulate a probability interpretation with a complex Boltzmann factor. 
However, there has been progress since several ways to work around the sign problem 
have been discovered, as will be discussed shortly. We propose a complementary idea 
which is to formulate the theory on a manifold with sufficiently small spatial volume 
such that perturbation theory is valid, and so that the calculation is valid at all 
temperatures and densities. Our particular choice of spatial manifold is S* 3 , which 
is motivated by the connection with the AdS/CFT correspondence. In a nut-shell, 
when the maximally super-symmetric M = 4 gauge theory is defined in a compact 
space, an S* 3 , its thermodynamics in the large N limit and at strong 't Hooft coupling 
can be addressed by the dual gravitational description. What is remarkable is that 
the phase structure seems to match on to the weak coupling description that is 
addressed in perturbation theory; namely, as the temperature is raised there is a 
confinement /deconfinement type transition, a phase transition that is described in 
the gravity dual as a Hawking-Page transition from AdS space to an AdS black hole 
space [1-3] . This kind of analysis has been extended to include an i?-charge chemical 
potential [4] . What we take from this is that it is interesting to investigate the phase 
structure of any gauge theory in finite volume, on an S* 3 , at weak coupling which 
is ensured if the size of the S 3 , R, is much smaller than the strong coupling scale 
R <C ^qcd- m this case a thermodynamic limit is ensured by taking the large N limit 
and genuine phase transitions occur. These transitions have all the characteristics of 
phase transitions that are expected in theories with finite N defined on flat space. 
The goal of the present paper is to extend this kind of analysis to SU(N) gauge 
theories with Nf fundamental quarks in the large N Veneziano limit [5] , i. e. with the 
ratio -jf fixed, with both finite temperature and baryon chemical potential. It will 
be interesting to relate our weak coupling results to strong coupling analysis based 
on the AdS/CFT correspondence. This will involve adding "flavour branes" to the 
basic set up to describe the quarks. 1 

1 The process of adding flavour to the basic AdS/CFT set-up has a huge literature. Most of this 
work addresses the case where the boundary theory is defined on R 3 and with fixed Nf, so that 
the flavours can be introduced in the probe approximation; however, the papers [6, 7] consider the 
case of global AdS with an S 3 boundary. 
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Our approach should have implications for understanding one of the outstand- 
ing problems in theoretical physics; namely, the behaviour of cold dense baryonic 
matter, which in essence corresponds to an understanding of QCD with non-zero 
baryon chemical potential //. Its resolution would permit contact between particle 
and nuclear physics via a quantitative description of bulk nuclear matter from first 
principles, and would set the study of compact stars on a firm theoretical footing, 
via the input of the QCD equation of state (energy density e(/i), pressure P(aO) into 
the Tolman-Oppenheimer-Volkoff equations for relativistic stellar structure. Indeed, 
such a programme is a necessary prerequisite for determining whether postulated 
ground states such as color-superconducting quark matter could ever exist in our 
universe. 

The current consensus [8] is that as baryon density n increases, the state best 
described as "nuclear matter", viz. a degenerate system of neutrons and protons 
with n ~ 0.45fm~ 3 which is the favoured ground state once \i exceeds its onset value 
\i — 924MeV, is somehow succeeded by an alternative degenerate system called 
"quark matter". The properties of quark matter have been the subject of intense 
speculation over the past decade; it has been suggested that for sufficiently low 
temperatures T, as a result of quark Cooper pair condensation at the Fermi sur- 
face, global, local and even translational symmetries may be spontaneously broken 
in the ground state, leading to exotic phenomena such as color superconductivity 
(CSC) or even crystallization [9]. From a theoretical standpoint these scenarios are 
most readily studied at weak coupling, which in the thermodynamic limit can only 
be quantitatively accurate at asymptotically high densities where ji ^> Aq C d, and 
yet the passage from nucleons to quarks referred to above clearly requires a non- 
perturbative treatment. At the very least, a reliable matching between perturbative 
and non-perturbative regimes is required. A recent calculation to this end has ap- 
peared in [10]. Other theoretical issues which naturally arise in this context include: 
what precisely is meant by "degenerate matter"? (this is usually taken to mean a 
system with a well-defined Fermi surface characterized by a momentum scale &f, but 
this definition is not gauge-invariant), and, does chiral symmetry restoration and/or 
deconfinement occur as /i increases? And, if so, to what extent do the transitions 
resemble those known to occur for /i = as T is raised? 

An important reason why a non-perturbative understanding of QCD with \x ^ 
has not progressed as much as in other areas is the unavailability of lattice gauge 
theory simulations performed using standard techniques [11]. A system with /i ^ 
is not invariant under time reversal, since a bias is introduced in favour of par- 
ticle propagation in the positive t direction. In the Euclidean metric this results 
in asymmetry between i and —i, so that for instance the Polyakov line defined by 
P(x) = TrU^Uo&t) has the property (P) ^ (P*)*, implying that in a medium 
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with n > the free energy of a static color source differs from that of an anti- 
source. Crucially, it implies that in general the Euclidean action is complex-valued, 
and hence has a fluctuating phase <fi = ^sS. Since the functional measure e~~ s is no 
longer positive definite, the Monte Carlo importance sampling used in lattice simu- 
lations is inoperable. Indeed, it appears that physically acceptable results can only 
be obtained if delicate cancellations are correctly handled over a much larger region 
of configuration space than that normally considered [12,13]. 

Attempts to evade this so-called Sign Problem fall into three classes. First, one 
can run simulations at /i — and attempt to calculate operator expectation values 
via re-weighting 

(O) = (1 1) 

where ((. . .)) denotes averaging with respect to a suitably chosen real measure. This 
approach has been found to be particularly effective in the vicinity of the quark- 
hadron phase transition for ^ ^ 1 [14], but must fail in the thermodynamic limit since 
the ratio of two partition functions ((e 1 ^)) ~ e~ cv , resulting in a disastrously poor 
signal to noise ratio as V — > oo. Other methods which can be applied in this physical 
regime, relevant for the hot medium produced in RHIC collisions, rely on analytic 
continuation of results generated from simulations with a real action, either at \x = 
by calculating successive terms in a Taylor expansion [15], or by simulating with 
imaginary chemical potential (corresponding to real constant abelian electrostatic 
potential) [16]. These latter approaches work in the thermodynamic limit, but are 
necessarily limited by a finite radius of convergence. The three methods have achieved 
some consensus in mapping out the quark-hadron transition line and determining the 
equation of state for small ^; whether they will ultimately prove capable of saying 
something about a possible critical point in the (//, T) plane is as yet unresolved. 

A second approach is to study gauge theories without a sign problem, ie. where 
the functional measure remains positive definite for /i ^ 0. These include QCD with 
isospin chemical potential (ie. with /ij = /i d — —fi u ) [17], and theories with real 
matter representations such as the fundamental of SU(2) (or the 6 of SU(A)) with 
Nf even, or any theory with adjoint quarks [18]. The generic feature of such models 
is a degeneracy between mesons and baryons at // = 0: hadron multiplets contain 
both mesonic qq states and qq or qq states which carry a non-zero baryon charge (in 
the case of fii ^ the latter role is played by 7r- mesons with non-zero In such 
theories the onset takes place for fi ~ M n , where M n is the mass of the pseudo- 
Goldstone boson associated with chiral symmetry breaking at /i — T — 0. Theories 
with sufficiently small quark mass so there is a large separation between Goldstone 
and hadronic scales are then well-described by chiral perturbation theory (^PT), in 
which only Goldstone degrees of freedom, some baryonic, are retained. For /x ^ /i 
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the resulting system is an (arbitrarily) dilute Bose-Einstein condensate (BEC) formed 
from tightly-bound qq bosons [18]. These models therefore fail to describe nuclear 
matter, but may still have something important to tell us about quark matter at 
higher densities. As /x is increased beyond onset there is a smooth rotation from the 
chiral condensate characteristic of the vacuum at T = /i = to a gauge invariant 
diquark condensate (qq) ^ 0. Since this condensate breaks a global, rather than a 
local, symmetry, the ground state is superfluid but is not a color superconductor. 
The predictions of for /i ^ /i D have been quantitatively confirmed by several 
lattice simulations [19]; more recent lattice works have explored the breakdown of 
%PT at larger values /i ~ ^qcd [20]. A related study of p-meson condensation in 
large- iV c QCD with p, ^ uses the methods of AdS/CFT duality [21]. 

We should mention two more radical approaches to simulating theories with 
p 7^ 0. It is possible to mitigate or even eliminate a sign problem by transforming to 
a different set of field variables via, eg. the exploitation of a duality symmetry [11]. 
Two recent papers have applied this idea to convert bosonic field theories in three 
dimensions to loop gases, which are then simulable [22,23]. A related treatment of 
nuclear matter in the strong-coupling limit of lattice QCD has also appeared in [24] . 
Finally, there is the possibility of generating representative field configurations by 
integrating a stochastic differential equation, the Langevin equation, in which the 
complex drift terms resulting from the action force the field variables to evolve in 
an extended, complexified space in which the large regions where observables are 
swamped by phase fluctuations can be avoided [25]. Again, this method has been 
successfully applied in certain bosonic cases [26]. 

In our approach, the fact that the action is complex for p ^ turns out not 
to be a problem but it does have important consequences. In the small volume 
theory, the effective action is defined over the eigenvalues of the the Polyakov line 
P, which is a unitary matrix whose eigenvalues can be written e t9i , % = 1,...,N. 
The functional integral reduces to an integral over the angles {^}. At large N the 
functional integral is dominated by a single saddle-point but since the action is not 
real this saddle-point configuration lies out in the complex plane where the 9i are no 
longer real. As a consequence (P) ^ (P 1 ")*. 

In the following sections we will summarize the derivation of the action of QCD 
formulated on S* 1 x S* 3 from one loop perturbation theory. Then we will present 
calculations of several observables as a function of the chemical potential, working 
at low temperatures, in both the small N, and the large N limits. For sufficiently 
small iV it is possible to calculate the partition function, and thus any observable 
derivable from it, by simply numerically performing the integrals over the gauge 
fields. Of course for finite N in finite volume there are no sharp phase transitions; 
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nevertheless, they show up qualitatively in the behaviour of observables. We present 
results for N = 3 which involves performing integrals over the 2 eigenvalues of the 
Polyakov line which is an SU(3) matrix. In this exploratory study we consider 
Nf = 1 Dirac fermion flavor, and consider the limit of a light quark mR = 0, and 
a heavy quark. What we find is that the fermion number, pressure, and energy rise 
in discrete levels as a function of the chemical potential. This is reminiscent of the 
quantum hall effect of QED [27, 28] where in our case the discrete levels result due 
to restricting our calculation to small spatial volumes. The magnetic field B which 
causes the quantum hall effect in QED is loosely analogous to our Rr 2 . The analogy 
can be taken a bit further in that in both QED with non-zero external magnetic field, 
and QCD in a small spatial volume both exhibit the level structure of the fermion 
number as a function of the chemical potential only in the low temperature limit. 
Increasing the temperature causes the levels to become smoothed out. Or, when 
B ~ Rr 2 is small (large volumes) the levels are also smoothed out. 

In the large N limit it is necessary to consider the gauge field, corresponding to 
the angles of the Polyakov line, as a distribution on a contour. From the equation-of- 
motion the saddle-point distribution of the Polyakov line eigenvalues can be calcu- 
lated analytically and plotted by mapping the angles from an arc on the unit circle 
to a contour over the same range of angles in the complex plane. What we observe 
is in agreement with the finite N results. The contour on which the Polyakov line 
eigenvalues are distributed is closed, corresponding to the confined phase, in between 
level transitions, and opens up while the transitions between levels are taking place. 
This is the characteristic feature of a third order, Gross- Witten transition [in the 
Ehrenfest classification], and indeed the third derivative of the grand potential is 
discontinuous at each level crossing [29-31]. 

2. Background 

Here we will summarize the one-loop formulation of QCD on S* 1 x S* 3 which was 
derived for SU(N) gauge theories with more general matter content in the beautiful 
paper [3]. This section does not present new material and is merely included for 
completeness. The partition function of QCD at finite temperature T = 1//3, for Nf 
quark flavours, each with a mass mj and coupled to a chemical potential ///is given 
in Euclidean space by 




(2.1) 



where ip and ip are the fundamental and anti-fundamental fermion fields, respectively, 
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and A is the SU(N) gauge field, = A^T a . The Lagrangian is 

Nf 

^qcd = ^Tr F (F, U F, U ) +J2$ f {Ip F {A) ~ loHf + m f ) ^, (2.2) 
with covariant derivative 

D^A) = M - A M , (2.3) 

and field tensor 

= [D,(A), D V {A)\ = duA, - d»A v + [A„, A v \ . (2.4) 

When the spacetime geometry is S 1 x S 3 , only the component Aq along the S* 1 has 
a zero mode. The idea is to construct a Wilsonian effective action for this mode. 
To this end we decompose Aq = a + gM) where, without-loss-of-generality, we can 
choose the background field a to consist of the diagonal elements of A while the 
fluctuation g£/ consists of the off-diagonal elements. The background field breaks 
the gauge symmetry from SU(N) to its maximal abelian group [/(l)^" 1 and gives 
mass to the off-diagonal modes which can then be integrated out. Gauge fixing with 
Feynman gauge and retaining the one-loop contributions puts the Lagrangian in the 
form 

^Qcn = - \<{D 2 Q {a) + AW)< " l^(Dl(a) + A^)A1 

- c(D 2 (a) + A«)c + $W F (a) - lo M + M)ip, 

where (ip) is an Nf component vector containing the ipf (?/>/), M (Ai) is an Nf x iV/ 
diagonal matrix containing the rrif (fif) as the diagonal elements, c and c are complex 
Grassmann- valued ghost fields resulting from gauge fixing. Here A^ and A^ repre- 
sent the scalar and vector Laplacians, respectively, where A^ = g^ 1 ^ 2 d^{g 1 ^ 2 d tl ) and 
A^A* = V j V j A i - R)A> with the Ricci tensor of S 3 . It is useful to decompose 
the spatial gauge field as Ai = Bi + C i: where Bi is the transverse (T) component 
with Vj-Bj = 0, and C» is the longitudinal (L) component with Cj = Vj/. Then we 
have 

^QCD = - \<{D 2 Q {a) + A«)< - is?(Dg(a) + A^)S? 

- ^C?(£>g(a) + A^)Cf - c(L» 2 (a) + A«)c + - 7o-M + Mty. 

(2.6) 

Performing the Gaussian integrals we obtain a simple form for the effective partition 
function (there is yet the integral over a to perform to obtain the full partition 
function), 

Z(a) = detail (-D 2 (a) - A< a >) det" 1 (-D 2 (a) - A^) det (fl F (a) - lo M + M) , 

(2.7) 
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where we note that on S 3 , the eigenfunctions of the scalar Laplacian have energies 
Eg and degeneracies df^ given by 

A^Y e (Q) = -s { e s)2 Y e (Q) , 

ef )2 = £(£ + 2)R~\ (2.8) 
4 s) = (£ + l) 2 , 

where £ = 0, 1, 2 . . . and R is the radius of S 3 . To obtain (|2.7|) we used the fact that 
the vector Laplacian acts on the longitudinal vectors Cj as 

A (u) (V fe /) = ( V 1 Vi5i - R{) V,/ = V fc ( V* VJ) (2.9) 

which results in the same spectrum as for scalars, A^ v ' L 'Ci = A^Cj, with the excep- 
tion that £ > 1 for vector fields. 2 This leads to the almost-cancellation between the 
M), Ci, c and c terms, with the £ = contribution, det^Q (— -Dq(oj) — A^), as the 
only remaining piece. 

To evaluate the fermion contribution we need the identity 

TV'ViVj = g ij V^j - -K. (2.10) 

where ft is the scalar curvature of S 3 . Then evaluation of the fermion determinant 
proceeds as follows performing the determinants over the flavor then spinor degrees 
of freedom, 



log Z/(a) = logdet(# F (a) - 7o-M + M) 
1 

~ 2 



logdet [- (# F (a) - 7oA*) +m J 



2\ f log (let - (D (a) - A*) 2 - A (/) + -ft + m 2 
The eigenvalues and degeneracies of the spinor Laplacian on S 3 are given by 

(a(/) - ift) v = -4^ , 

ef 2 =(£ + \) 2 R~\ ( 2 ' 12 ) 
dtp =£(£+!), 

where £ — 1, 2, . . .. Regarding the S 11 around which the fermions have anti-periodic 
(thermal) boundary conditions, the eigenvalues of the operator D (a) are discretized 
in terms of the Matsubara frequencies, u n = (2n + 1)tt/(3, and given by 

D (a) ^iu n -a a T a . (2.13) 



2 Thc derivative of the t — mode vanishes since it is a constant. 
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The fermion contribution then takes the form 



log Z,(a) = 2N f Tr R £ £ [log L* n + {ef m) af) 

neze=i L v ' (2.14) 

+iog(^ t + (4 / ' m) + /^ + «) 2 ) 



where e^/'" 1 ^ = y e^p 2 + m 2 . 

We define the Polyakov line order parameter, P, by the path-ordered exponential 
of the temporal gauge field. In terms of the constant temporal background field 
a = i0/ (3, it is 

P = & e lg*M*) = e P<* = e M = diag(e*\ ...,e idN ) . (2.15) 
Then, following [32], 

oo oo / 1 s n N 

1=1 n =l 71 i=l 

(2.16) 

where we have dropped terms that are independent of the gauge field as their con- 
tribution cancels when taking expectation values. 

We follow a similar procedure to obtain the boson contribution. Given the 
eigenvalues and degeneracies of the transverse vector Laplacian on S 3 , 

A^B}(£l) = -ef' T)2 B}(Cl) , 

e^ 2 = (£ + l) 2 R- 2 , (2.17) 
dp T) = 2£(i + 2) , 

for £ = 1, 2, . . ., we have 

log Z b (6i) = ^logdet, =0 {-D 2 (a) - A«) - logdet (-D*(a) - A^) 

oo / oo \ N (2.18) 

n=l \ 1=1 J ij=l 

where we used the fact that the trace in the adjoint representation 

N 

Tr A (P) = cos(n(6i - 9,)) . (2.19) 

ij=l 
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Adding the boson and fermion contribution the total one-loop effective action is 



Sfa) = jr - (1 - z b (nf3/R)) cosing - 6,)) 

n=l i,j=l 

oo / -.w N 

+ ~ -N f z f {nP/R, mR) ^ [ e n ^ +in ^ + e ~^-^] 



n=l 

where we have defined 



(2.20) 



z b {P/R) = Y,d ( t T) e-^ T) 

oo 



2^£(£ + 2)e-^ +1 ^ (2.21) 



and 



z f {(3/R,mR) = J2 d e f)e 



=i 

oo 



i=i (2.22) 

2m 2 R 3 mR 

-K 2 (pm) —Ki(pm) 



P yr ' 2 

r-oo ™2 I 1 



✓•OO .4 _|_ J_ 

+ 4 / da; 2™ 4 -, sin(/3Vx 2 - m 2 R 2 /R) . 

JmR e nx + 1 



Notice that the first term in the sum ( [2. 20 ), up to an unimportant constant 



can also be interpreted as the Vandermonde Jacobian contribution resulting when 
converting from an integral over unitary matrices of the form P = e ld to an integral 
over the eigenvalues, of 0. That is, 



N N /an 
sin 2 ' 3 



dp= I n^n 

i=l j<i 
N 

Vdm 

i=l 



(2.23) 

j exp [ — S\ 



In the low temperature limit (/3 — > oo) we have z&(oo) = and so the gluonic 
contribution to the action reduces to the Vandermonde piece, SVdm, defined above: 

S = SVdm + S f (2.24) 



- 9 - 



and so in this limit the theory reduces to that to that of an N x N unitary matrix 
model 

Z = JdPe- NTrV ^ , (2.25) 
with a potential determined by the quark contributions that we can write as 

oo 

V{P) = -^2<T t [log(l + e^- £e) P) + log(l + efl-^pt)] . (2.26) 

£=1 

with the definitions that we use from now on 



a e = 2£(i + 1)^ , e e = ef m) = j m * + (I + \y R -2 . (2.27) 



3. Finite N: N = 3 



Even though there are no sharply-defined phase transitions for a finite iV theory valid 
in a small spatial volume, taking the low temperature (large S 1 ) limit is not that 
far removed from a limit where well-defined phase transitions are possible: as the 
temperature is decreased transitions as a function of the chemical potential appear 
more and more sharp, even if on a microscopic scale they are always continuous. To 
develop an understanding of the physics at non-zero chemical potential we calculate 
several observables from the low temperature partition function, 



r r 00 1 

Z(f3/R)= [d6}exp - J2~ [Tr A (P") + (-l) n N f z f (n/3/R,mR)e n ^Ti F (P n )] 

J [ n=l 71 

(3.1) 

where [d9] = YliLi d&i- These include: 
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1 f d lo Z \ 

Fermion number: ^ = — I — ) (3.2) 

P \ dft J 

Polyakov lines: ^ = ^ J [d9] e~ s ^ e i0i j (3.3) 

^-i = | / [^]e- 5 (^e-^) (3.4) 

Energy: £ = -7^ + //^ (3.6) 

Chiral condensate: (V'V') = — tttt nm I — I (3.7) 

p V3 m^o y dm J 

Average phase: (e -) P q = (3.8) 



where the pq in the average phase refers to the phase-quenched theory to be discussed 
later. In what follows we first calculate each of these observables for arbitrary N, not 
performing the integrals over 0j. Then we present results for N = 3 by integrating 
numerically over the 8{. Each observable is calculated as an expectation value with 
the form 

= / [dfl] e- s O fM^e-SQ 

J[d6]e- S ^3 JdO^e- 3 K ' 

where the integrals over 8 3 drop out as 9 3 = —Q\ — 9 2 by the SU(N) condition. 
These results are multiplied by factors of (3 and / or V3 as needed to make them 
dimensionless. In this paper we present results for Nf = 1 Dirac fermion flavor. We 
calculate the above observables considering first the case of a massless quark, and 
then the case of a quark with large mass. 



3.0.1 Fermion number jV for m = 



The fermion number jV gives the number of quarks minus the number of antiquarks 
in the volume of S 3 , V3 = 2tt 2 R 3 , where R is the radius of S 3 . From Figure |], which 
shows jY as a function of [iR for low temperatures for a single massless quark flavor, 
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Figure 1: Expectation value of the fermion number as a function of the quark chemical 
potential for QCD on 5 1 x S 3 . TV = 3, N f = 1, m = 0, p/R = 30 (low T). 



an occupation level structure is apparent. The fermion number in this limit is 



1 /<91ogZ 



(3.10) 



^f/i^- s £X> ( £+i) 

7 £=1 i=l 



e Pn + e -i0i+/3(^+i)/R 



where the derivative of the action with respect to the chemical potential brings down 
a factor of nj3, leading to a geometric series which gets summed to give the Fermi- 
Dirac distribution. Ignoring the Polyakov line for the moment, it is clear from the 
general form of the Fermi-Dirac distribution function, 

f(e e ) = - \. r, (3.11) 

that the transitions occur when ee — \i changes sign, i.e., when /i passes an energy 
level. When /j, <C E\ then / ~ 0. When /i ^> si then / ~ 1. The system is in 
transition for \/3(ee — small. This shows that each level L has a net number of 
quarks given by 



jV l = NN 



L 

^2^+1) 



(3.12) 



- 12 - 



Each new level starts at 



(fiR) = L + ±, (3.13) 

has a level width 

(A^) A i = 1 , (3.14) 

and a transition width given by the width of the Fermi-Dirac distribution function. 
This goes like the temperature, 

(AnR) L ~ ^. (3.15) 



Since the level width in fx goes like 1/R, then the levels should become narrower 
with increasing spatial volume when considered as a function of the chemical potential 
alone. In addition, as we will show in the next section, the width of the levels also 
depends on the quark mass m, and taking m large also causes the steps to become 
narrower, at least until /i >> m. For large values of fi or small values of f3 (finite 
temperature) the transition width is larger. 

It is interesting to compare with the results of [23]. The authors observe the 
same level structure in the particle number in the non-linear 0(2) sigma model. 
Their results for the particle number as a function of the chemical potential (not 
multiplied by the spatial extent) show that the levels become narrower as the spatial 
volume is increased, and the particle number appears more continuous. Our results 
indicate a level width A/z = ^ (in the massless limit) and are thus qualitatively 
consistent with theirs. 



3.0.2 Polyakov lines: and for m = 



At zero chemical potential, = (TrP) and = (TrP" 1 ") are complex conju- 

gates, but for non-zero chemical potential this is no longer true 3 . The Polyakov line 
expectation values are 

1 r N 

&>i = ^ [d9]e-' s ^e ie % (3.16) 
J i=l 

and 

N 

= -= / [d6]e- s J2 e ~ i6i - ( 3 - 17 ) 

J i=l 

3 The matrix model studied in [33] has a similar fermion term and also shows this effect. 
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Figure 2: Expectation values of Polyakov lines &\ and as a function of the chemical 
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Figure 3: and .^Li as a function of the chemical potential at the first transition 
(Left), and the fourth (Right). JV = 3, N f = 1, m = 0, (5/R = 30 (low T). The width of 
the deconfined regions increases with fiR. 



Figure || shows &\ and as a function of /ii?. Each spike in and 

corresponds to a level transition in jY . Even though their behaviour as a function 
of fiR is similar, the peaks of always preceed those of at the start and finish 
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12345 12345 

Figure 4: (Left) Pressure and (Right) energy for N = 3, Nf = 1, m = 0, (3/R = 30 (low 
T). 

of each level transition. 

Taking the temperature lower than that of Figure |2] causes the peaks to be 
narrower while maintaining the same height. Taking T — > would make the peaks 
appear as infinitely narrow lines, occurring precisely at fiR = 1.5, 2.5, .... 

In Figure [| we compare the first and fourth transitions. For non-zero T, as 
fiR increases the transition width, and thus the width of the deconfined regions, 
increases. 



3.0.3 Pressure V for m — 



The pressure indicates how the system responds to changes in the spatial volume. 
The pressure multiplied by the 4-volume is shown in Figure § (Left) for N = 3. In 
the p — > oo, m — > limit the expectation value of the pressure is given by 



_ 1 ( dhgZ 
P V dV 3 



-R 

3(3V 3 Z 



/3->co 3RV 3 Z 



„ oo N 

J [d^]e- s ^^2^+ !)(£+!) 



(3.18) 



_ e #i _|_ e 



-i0i+p(t+i)/R 
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which indicates that each level has a pressure 

"3 1 



(3.19) 



3.0.4 Energy E for m = 

The energy can be evaluated from the pressure and fermion number above and is 
plotted in Figure £| (Right). The energy is calculated as 

E = -VV 3 + h<jV 



/3— >oo 

N f 



oo N 



o i „■ — 1 



e=i i=i 



(3.20) 



e /?M _|_ e -iBi+l3(e+±)/R 



This shows that the energy levels are not horizontal. The factor of /i in front of the 
fermion number in the first line causes the levels to rise linearly with /i. The energy 
of each level is given by 



E L = NNf W + 1) (A* - W + \)I R ) ■ 



(3.21) 



i=i 



3.0.5 Normalized , V, E for m = 0; the Stefan-Boltzmann limit 



The Stefan-Boltzmann limit is the zero interaction free fermion limit. On S 1 x S 3 
we obtain it from the one- loop result taking all the 0j = 0, corresponding to the 
deconfined phase, e.g., for the fermion number 



SB 



NN f 2£(£ - 



e=i 



M+h/R 



(3.22) 



The normalized fermion number, jV j ' jYsb, is shown in Figure |5|. The behavior of 
V jVsB an d EjE SB as a function of /i is almost indistinguishable from -jV j ,JV SB . 
These results might at first seem at odds with the results for and <^_i, in the 
confined regions, since the Oi are set to in Vsb, ^sb, and Esb- The resolution is 
that sufficiently far within the confined regions the observables are independent of 
the 9i. When fiR is not close to £ + | for / = 1,2,..., then in the /3 — > oo limit the 0, 
terms drop out of the exponentials in the Fermi-Dirac distributions, 



(3.23) 
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Figure 5: Average fermion number normalized by its Stefan-Boltzmann value as a function 
of the chemical potential. N = 3, m = 0, fi/R = 30 (low T). 
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Figure 6: Average phase as a function of chemical potential. = 3, m = 0, j3/R = 30 
(low T). The average phase only differs from 1 when &\ / &*_ v 

3.0.6 Average Phase (e l<t> ) pq for m = 



The average phase is important in that it shows where the sign problem is severe. It 
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is given by 



(e^)p 3 - (e Jlm(5) ) pg - r^] e R e (-5)' (3.24) 

where the denominator is the "phase quenched" (real action) partition function, 

Zpq = J [d6] |e~ s | = J [d^]e Re( - s) . (3.25) 

Numerical results for the average phase are presented in Figure ||. The results suggest 
that the sign problem increases in severity with the chemical potential, then levels 
off (large fiR results suggest it levels off near 0.83). However, the average phase is 
only different from 1 while &\ and gP*_ x differ. This explains why the average phase 
is always 1 in the middle of a transition. In fact, We observe that {e l ^) pq is smallest 
(largest) when \ — is largest (smallest). 

3.1 Chiral condensate (?/n/>) 

The chiral condensate is given by 
(ipip) = — — — lim ^ 



/3Vs m-^o \ dm 
1 

w?. 



N f m /•,,„, a^J^£(£+l) 



->■ 



^oo.m^o n 2 R 2 J ^ ^ (£ + i) 

i=\ i=l v 



e /3v _|_ e -t6i+P(i+±)/R 



In the massless limit (ipip) = as expected since it is a perturbative result, but it is 
interesting to note that in the light mass limit is linear in m. 

3.2 Continuum results (mR — > oo) 

Since all of our observables are functions of the dimensionless quantities (3/R, mR, 
or /j,R, then we can obtain a "continuum" limit by taking one of the following: 

• (3/R small (high temperature perturbation theory), 

• \iR large (high density perturbation theory), 

• mR large (heavy quarks). 
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We take mR large to maintain validity in the region of the confinement-deconfinement 
transitions. For sufficiently large mass the fermion contribution to the action can be 
converted from a sum to an integral using the Abel-Plana formula and then simplified 
by taking mR — > oo. This proceeds as 

oo 

z f (np/R, mR) = 2j2 W + l) e - n V™ 2 +('+§) 2 *- 2 

I'OO 

= 2 dy (y> - I) e-tV^* 
Jo 

+ 4 r d yi^ sin ( n P vV - m2R2 / R ) 

JmR e f + i \ / 



/■oo 

> 2 / dy (y 2 - i) e -nPVy 2 +™ 2 RVR. 

miJ->oo Jq 

The resulting integral formula removes all the periodic structure from the observables. 
It should be noted that this result is only valid in the vicinity of, fxR = mR. Once 
the chemical potential becomes sufficiently large compared to the fermion mass then 
the periodic structure returns and the system behaves again as it did in the massless 
case. 



3.2.1 JV for m ->■ oo 



For QCD with non-zero quark mass the expectation values <jV , V, and E exhibit 
"Silver Blaze" behavior: bulk observables are zero until onset [12] which occurs when 
the chemical potential reaches the value of the lightest quark mass. Here onset occurs 
for ji — m. The fermion number is given by 



mR->-oo IV f 



/3-s>oo 



-j- J me- S J o dy2(y 2 -l/4)g 



(3.28) 



Each level L has height 



h L = 2NN f J2 i (?+ l ) 



h y = 2NN f J dy(y 2 -l), (3.29) 



and level width 

(AfiR) Al = ^(L + i + l) 2 + m 2 R 2 - \/(L + \) 2 + m 2 R 2 
->■ \JJy + dy) 2 + m 2 R 2 - \/y 2 + m 2 i? 2 



(3.30) 



->■ 0. 
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Figure 7: Fermion number: (Left) unnormalized, and (Right) normalized by its Stefan- 
Boltzmann value as a function of the chemical potential for large quark mass near onset at 
/iR = mR = 30. N = 3, Nf = 1, j3/R = 30 (low T). The dots are calculated using the full 
sum form of Zf. The curves are from the integral approximation. 



The Stefan-Boltzmann values shift accordingly as well, e.g., for the fermion number 

P 0v 



SB 



NN f J2 2 Kl + l ) 



j3— >oo,mR— >oo 



-> NN 



dy2(y< 



T?R 2 



(3.31) 



The fermion number for large quark mass is plotted in Figure 0. The dots are 
calculated using the full sum form of Zf (line 1 of eq. ( |3.27| )) and the green curves 
are from the integral approximation (line 4). As expected jV is close to until 
the onset transition at fi — m. The smoothness of jV as it increases from onset 
results from the decreased level width in the large mR limit. It is only temporary 
though; as fiR is increased further from mR the levels spread out and eventually the 
observables behave as they would for mR = 0. This is reflected in the oscillations 
that develop as fiR is increased when considering the full sum form of Zf. The lack of 
these oscillations in the integral form shows that this approximation breaks down for 
values of fiR away from the onset transition at mR. How long this breakdown takes 
depends on mR. The larger the mR, the further in fiR we can go before breakdown. 



3.2.2 Polyakov lines: £P\ and <^_i and average phase (e l ^) pq for m — > oo 



As in the mR = discrete case, in the large mR limit the behavior of <^_i precedes 
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Figure 8: (Left) &\ and as a function of chemical potential for large quark mass 
near onset at /iR = mR = 30. N = 3, Nj = 1 (3/R = 30 (low T). The dots are calculated 
using the full sum form of zj. The curves are from the integral approximation. (Right) 
The average phase is calculated under the same conditions. 



that of 0P\ as a function of fiR. The transition in fiR occurs around onset at mR 
and appears sharper for larger mR. Results for mR = 30 near the transition are 
presented in Figure [8] (Left), which shows and 0P-\ using the sum form of Zf 
(black dots) and the integral approximation (green and blue curves). It is clear again 
that the integral approximation breaks down after the transition. 

Figure H (Right) shows the average phase near the transition. After the tran- 
sition, and before fiR is sufficiently large that oscillations return full force, there is 
a brief respite from the sign problem in that the average phase is close to 1. The 
small spike in {e l ^) pg at /i = m corresponds to the value of fiR where &\ and 1 
differ maximally and serves as a good indicator of the location of the transition. 
Eventually, regardless of how large we make mR, the oscillations always return by 
taking sufficiently bigger. 



4. The Large TV Theory at Low T 



In this section we turn to an analysis of the large- iV theory, specifically in the low- 
temperature limit T R" 1 . We showed in ( |2.25| ) that in this limit the system 
reduced to a unitary matrix model with a potential determined by the quarks. 4 

4 A similar matrix model but defined in the canonical ensemble and truncated so as to include 
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Without the chemical potential, the potential V(P) in ( 2.25 ) would vanish in 
the low temperature limit. In this case, the Vandermonde piece dominates leading 
to a repulsion of the eigenvalues and so the integral is dominated by a uniform 
distribution of the angles 0j around the circle. Consequently (TrP n ) = and the 
theory is in the confining phase. However, if as T — > we simultaneously tune 
//—>■£ keeping /3(n — e) fixed, where e is one of the energy levels of the fermion 
system then the fermion term can compete with the van de Monde piece and a phase 
transition can occur. In fact this is a version of the Gross- Witten transition [29] (see 
also [30,31]). The Gross- Witten transition is typically a third order phase transition 
that occurs in large iV gauge theories in finite volume and involves a transition from 
a configuration where the density of eigenvalues of the Polyakov line {e l9i } lies on 
the unit circle in the complex z = e td plane to a configuration with a gap. We 
loosely refer to this as a confinement /deconfinement transition because the confined 
phase defined by a distribution that lies on a closed contour is smoothly connected 
to the phase where the eigenvalues are uniformly distributed around the unit circle 
and (TrP n ) = 0, for n / 0, indicative that it costs an infinite amount of energy 
to propagate colour charges. Correspondingly the deconfined phase defined by a 
distribution that has a gap is smoothly connected to the configurations where all 
the eigenvalues e %di = 1, i.e. (TrP n ) = N where colour charges are free to propagate 
indicative of a plasma phase. We will find that, as suggested by the N = 3 results, 
increasing the chemical potential induces a series of such transitions, the novel aspect 
being that the eigenvalues no longer lie on the unit circle due to the complex action. 

Without-loss-of-generality we will choose /i to be positive so that only the con- 
tribution from the quarks survive and the quark potential takes the form 

V(P) = - log (l + e^-^P) . (41) 

i 

Notice that the potential is not hermitian. This is the usual "sign problem" in the 
presence of a chemical potential. In the context of the matrix model, it means that 
the saddle point that dominates in the large N limit will lie at complex angles $i. 
In this respect the matrix model is of the "holomorphic" type that appears in the 
Dijkgraaf-Vafa approach to supersymmetric gauge theories and in particular [35] 
which considered the unitary version. 



4.1 Single level model 

For the moment let us focus on what can happen with a single level e = S\ and 

a = <j\. Defining an effective fugacity £ = e^~ e > the effective action on the angles, 

only the n = 1 term was considered in [34] . We work directly in the grand canonical ensemble and 
it is important in our analysis that we do not truncate the sum over n in the fermionic sector. 
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including the Vandermonde piece, is 

1 N rO - \ N 

i,j=l i=l 

(4.2) 

where we have added the Lagrange multiplier M to enforce the detP = 1 constraint, 
i.e. Ei^ = 0. 

In the large N limit, the integral over the angles is dominated by a saddle point 
obtained by solving the equation-of-motion that follows from (|4.2|) 



1 + £ e ld > N ^ V 2 J 



(4.3) 



As we remarked above, our system is a unitary matrix model with a non-standard 
kind of potential term V(6) coming from the fermions. What is novel in the present 
context is that the potential is not real and as a consequence the saddle-point con- 
figuration will lie out in the complex plane. In fact if we define z-i = e Wi then in the 
presence of the non-real potential the Zi will move off the unit circle in the z-plane. 
As a consequence if we define the observables 

1 N 

0» n = (TrP n ) = -Y,e <nei > ( 4 - 4 ) 

we will find that 3^- n ^ &* n on the saddle point solution. 

Before we plunge in and solve the matrix model, it is useful to note that when £ 
is either very small or large, the potential vanishes (in the latter case following from 
using the constraint £\ 9i = 0) and so we expect the {z^} to be uniformly distributed 
around the unit circle. Taking the equation-of-motion in terms of the variables 



cr^Zi 1 Zi + Zj 



(4.5) 



it then follows that as £ — > the Lagrange multiplier M = while as £ — > oo 
the Lagrange multiplier M = a. The Lagrange multiplier M has an important 
interpretation following from 



. 1 \ at;Zi T d log Z 

M =N\Y^r^^r (4 - 6) 

and so M is the effective fermion number, jL. So as ji varies from ji <C e to ^> e the 
picture is that the energy level becomes occupied and the effective fermion number 
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jumps by the factor a. The question before us is to establish how this transition 
occurs. 



The small £ confined phase 

We make the hypothesis that as £ increases from that the eigenvalues are 
continuously distributed around a closed contour C in the z-plane at least up to some 
finite value of £. In the large N limit, we can describe the distribution of eigenvalues 
by an analytic function g(z) for which g{z)dz along C is real and positive. It is useful 
to think of g(z) in terms of a conformal map of the cylinder z(s), —tt < s < ir, s G R, 
for which the eigenvalues are uniformly distributed along the real axis between — 7r 
and 7r, to the z-plane: 

i A r ds r dz 



i=i 



where the contour is obtained from the inverse map z(s) obtained by solving the 
differential equation 

<§ = M , (4.8) 

subject to the initial condition that when £ = we have z = e ls . For consistency we 
must also have z{— tt) = z(ir) so that C is closed. The normalization condition is 

dz , . , 

m<m = i , («) 

where the left-hand side can be evaluated using Cauchy's theorem. The SU(N) 
condition, 9i = 0, translates into the constraint 

j^g( z )\ogz = , (4.10) 

where the branch cut of the logarithm is taken through the contour C at the point 
z(±tt). 

In the large- N limit the saddle-point equation (|4.5|) becomes for z G C 

= 0^1^1=7 • 2rw ^-lffe- (4 - n) 

The here indicates a principal value which is the required prescription given that 
z lies on integration contour C. In practical terms this means that the right-hand 
side is the average of two terms with z infinitesimally just inside and just outside the 
integration contour. The right-hand side can then be evaluated by Cauchy's theorem 
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yielding an expression of the form —zg{z) + • ■ ■ where the ellipsis represent other 
terms that arise when g(z) has other poles inside C. It is clear that g(z) can only 
have poles at z — — | and 0. The condition ( }4.9| ) implies that the pole at must be 
inside C and have unit residue. The only question is whether the pole at z = — | is 
inside C or not. In the small £ region, it will be outside — indeed this will define the 
small £ region. Consequently, from 



r dz_*Zz e {z) = 



The form of g(z) is then completely fixed by (|4.1 1|) to be 



One can check that the constraint ( f4.10| ) is satisfied. Notice that as £ — > 0, g(z) — > -, 



which is the uniform distribution around the unit circle. The contour C follows by 
solving (|4.8| ) for s to give 

e is = z{l + izy (4.14) 

and then inverting to give z(s). Notice that the integration constant is fixed by 
requiring that as £ — > we have z = e ls (the uniform distribution of eigenvalues). In 
this small £ confining phase the effective fermion number vanishes M = and the 
Polyakov line expectation values are 

As advertised, as a symptom of having a complex action, ^ 

The question is what happens to this phase as £ increases. The issue is that 
the pole of g{z) at z — — j must lie outside C, however, this cannot be maintained 
when £ becomes large enough. To see what happens notice that g(z) vanishes at 
z = ~ a P°i n t that lies outside C for small £. As £ increases to the value 

f = &=(T^s; < 4 - 16 > 

g(z) vanishes precisely on C at the point z(±tt) on the negative real axis. At this 



point C develops a kink and beyond this value of £ the inverse solution z(s) to ( [4.14j ) 
ceases to be closed z{— tt) ^ z(tt). This signals that a phase transition will occur 
to a configuration where the contour opens into an arc, just as in the matrix model 
solved by Gross and Witten [29] . The line of the phase transitions in the (/z, T) plane 
corresponds to the straight line 

fi = e-T[(l + a)log(l + a)-a\oga] , (4.17) 
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valid in the low temperature limit. 



The large £ confined phase 

We can solve for the large £ phase in a similar way. In this case the pole of q(z) 
at z = — | is now inside C. The solution turns out to be 



which satisfies ( 4 . 1 1| ) and ( |4.9| ) with J\f = a. In addition, one can check that the 
constraint ( |4.10| ) is satisfied. Notice that as £ — > oo, we have 0(2) — >■ -, as required. 



The contour C follows by solving (|4.8|) to get 



e is = ^— . (4.19) 

In this phase the effective fermion number M = a — the level is now occupied — and 
the Polyakov line expectation values are 

^i = |» ^-i = 0. (4.20) 



So comparing with ( |4.15 ), the behaviour of swaps over along with the replace- 
ment £ — > 

As in the small £ phase, the large £ phase persists until the zero of q(z), at 
z = which lies inside C for large enough £, just touches C. This occurs when 

f = 6= (i±^. (4.21) 
a' 7 

For smaller values of £, the contour C is not closed and the phase does not exist. 
Notice that the points of transition £ = £1 and £ = £2 satisfy £1^2 = 1- In the (ji,T) 
plane the boundary lies along the straight line 

H = e + T[(l+a) log(l + a)-a log a] , (4.22) 

valid in the low temperature limit. 

The deconhned phase 

In the region £1 < £ < £ 2 , experience with the Gross- Witten matrix model 
suggests that the eigenvalues lie on an open contour C. In this case, the model can 
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be solved by the standard resolvant/spectral curve method. The resolvant is defined 
as the function 



«>(*) = -tfY— 

v ' N z — Zj 
j J 



(4.23) 



By hypothesis, in the large N limit the eigenvalues distribute themselves on an open 
contour C and oj(z) is a function which is analytic everywhere in the z-plane except 
along a square- root branch cut running along C with branch points, z and z*, lying 
at the endpoints, as illustrated in Figure || The eigenvalues are then distributed 
along the cut between the branch points, such that 

co(z) = - [ ^(z>) Z -±^. (4.24) 



c 



it follows from this representation that as \z\ — > and oo 



lim u(z) = 1 , lim u(z) = -1 . (4.25) 

\z\— >0 \z\^oo y ' 



From the Plemelj formulae [36], the equation-of- motion ( 4.11|) is simply the con- 
dition 

zV'(z) = -~[cj(z + e) +u( z -e)] , zeC, (4.26) 

where e is an infinitesimal such that z ± e lies on either side of the cut. 5 It follows 
from ( f4.24Q and the Plemelj formulae that the density of eigenvalues is obtained as 
the discontinuity of u)(z) across the cut C: 

zq{z) = - [u{z + e) - cj(z - e)] , zeC . (4.27) 
In particular, this means that an average of the form 

can be used where C is a contour which encircles C as illustrated in Figure ^. This 
can then be evaluated by pulling the closed contour C off the cut to pick up the 
residues in the z-plane. 

The resolvant can be written down based on the following conditions: it is an 
analytic function in the z-plane apart from a single square- root branch cut, it satisfies 



3 This is a definition of the principal value in (4.11). 



-27- 



Figure 9: The resolvant uj(z) is naturally denned on the cut complex z-plane with the 
eigenvalues taking support along the cut C. 



the conditions (|4.25 ), and g(z) = g(z*) = 0. The solution to the equation-of-motion 
((126]) is of the form 



u (z) = -zV'(z) + f( z )y/( z -z)(z-z*) , zg(z) = f(z)y/(z-z)(z-z*) . 

(4.29) 



The endpoints z, z* and the function f(z) are then determined by imposing ( [4. 25| ) 
and requiring that u)(z) is regular at z — — \. These conditions uniquely determine 



/(*) 



(4.30) 



and the endpoints of C are given by 
-1 



M 2 + 1 + a - Afa + 2t v / M(a-Af) (1 + a) 



(4.31) 



What remains is to fix M by imposing the final SU(N) condition (|4.10|) which be- 
comes 



f dz ^ . . Yogz = . 
Jq Aniz 



(4.32) 



Pulling the contour off the cut, deforming it around the poles at z = and z — — j 
collecting their residues, and adding the discontinuity from the branch cut of log z 
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gives 



■ u(z) log z = lim 



■ u(z) log z — (f oj(z) log z 



o 2iriz 27T22; 

_ e (4.33) 

+ f —— co(z)logz + / — u(z) = 0. 



77 — yC/O 



2niz J_ v z 

where e and rj are cut-offs and the integrals around z = and oo are defined on the 
contours ee td and rje 10 , < 9 < 2tt, respectively. The limit is well defined since the 
factors of log e and log rj separately cancel. What remains is the condition 

c (^-A/y-M + A/y+^ 

5 M"(1 + <t-M) 1 +°- j * ' K 1 

This equation determines M as a function of £. An important check on the solution is 
how the deconfined phase interfaces with the small and large £ closed phases. When 
A/" = we have £ = £i and at this point ImS = and so the two branch points come 
together and one can readily verify that the density q(z) matches that of the small 
£ phase at the edge of its existence at £ = £i. The same happens at Af = a when 
£ = £2 where the gapped phase merges continuously with the large £ closed phase. 



From ( 4.34 ) it follows that across the transitions at £ = £1 and £ = £2, M and 



its first derivative are continuous, however higher derivatives are discontinuous. 
Since M is the effective fermion number, the first derivative of the grand potential, 
it follows that the transitions are third order just as in the original Gross- Witten 
model [29]. The behaviour of the effective fermion number as \x is increased past e 



is shown in Figure [TO . 



The expectation values of the Polyakov line in the deconfined phase can be 
determined from an expansion of the resolvant for z — > and z — > 00 [38], 

00 j 

u;(z) = -l-2^-^ n , (4.35) 

n=l Z 

00 

u{z) = l + 2^z n 0>- n . (4.36) 

n=l 

where n is the number of windings. For a single winding, the Polyakov lines are 

& x= M 1 & 1 = ^t (4 37) 

where A/" = A"(£) via the inversion of (|4.34[) and £ = £ \ The behaviour of &±x 



is shown in Figure O. Notice how £?±x vanish in the small and large £ confined 
phases, respectively. 

6 For another example see [37]. 
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Figure 10: The effective fermion number across the pair of Gross- Witten transitions from 
the small £ confined phase through the deconfined phase to the large £ confined phase. 




The behaviour of the contour C is shown in the three plots that make up Figure 
12| . These plots also show how the pole at z = — | is "eaten" as one goes from the 



small to large £ phases. Finally Figure shows schematically the phase diagram in 
the (n, T) plane in the vicinity of fi = e and at low temperature. 
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Figure 12: The contour C, which gives the distribution of the eigenvalues of the Polyakov 
line, showing the transition from the small £ closed phase (in red), the open phase (in blue) 
and the large £ closed phase (green). 

4.2 Multi-level model 

In the full QCD model there are, of course a series of energy levels as in (|4.1|) . For 
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Figure 13: Phase structure in the single level model. For finite T as n is increased 
there is a Gross- Witten transition to the deconfined phase and then another Gross- Witten 
transition back to the confining phase at a larger /i. In the process the pole at z = — | 
goes from outside the contour C to inside. 



a series of levels we can write 

V'(0)=iA/'0-5>*log(l + &e*) , 



with 



6 = e 



(4.38) 



(4.39) 



In general a matrix model with such a potential could exhibit a rich set of 
phases where the eigenvalues have support over a multiple set of open contours, 
the "multi-cut" solutions. These multi-cut solutions arise when the potential has 
multiple degenerate minima and the cuts are located around the minima. In the 
present case the potential apparently does not exhibit multiple minima and so we 
suspect that these multi-cut solutions never dominate the ensemble. Henceforth, 
we will only consider solutions with a closed contour, the confining phase, and with 
single cuts, the deconfined phase, leaving a detailed analysis of multi-cut solutions 
for the future. 



When T is low enough, experience with the single-level model suggests that as /i 
increases the system goes through a series of phase transitions with a confined phase 
lying in a window of /i with < jj, < E£ + i, for which the closed contour C contains 
the poles £ K , with k < t. In this case the saddle-point equation is solved by 

A/" + 1 , inOn 

(4.40) 



g(z) 



+ 
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T 




s ' S 2 M 

Figure 14: A schematic picture of the phase structure of the multi-level model. The 3 
portions of confining phase are characterized by the fact that the closed contour C contains 
the poles at {0}, {0, — ^-} and {0, — — respectively, as fi increases. 



with the effective fermion number 



K<£ 



(4.41) 



indicating the first I levels are filled. 

The actual distribution of eigenvalues follows from integrating ( |4.8|) which gives 

The phase only exists over the finite region of \x for which z(ir) = z(—ir) so that 
the contour C is actually closed. The boundary of this region can be obtained by 
the condition that on the boundary q(z) = at some point on z G C. This defines 
a triangular shaped region in the (fi, T) plane which at T = covers the region 
£i < A* — £ e+i- The apex of the triangle lies at a special point where g(z) has a 
double zero on C. So with a number of energy levels the boundary of the confined 
phases describes a saw-tooth pattern in the (/i, T) plane as illustrated in Figure TA 



4.3 The continuum model 



In the limit, m 3> R" 1 , the fermionic levels form an approximate continuum starting 
at e = m and extending up to roughly 2m where the discrete structure starts to 
manifest itself again. As long as the temperature is not too low, T ^> (along 
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with T <C required to have a unitarity matrix model) the sum over levels can be 
approximated by a continuum of the form 



17 /"OO 

7(0) =iM9--= dy ^log (1 + £e-V e ) 
v 71 " </o 



(4.43) 



where £ = e 13 ^ m ^ and a = y/2TTjf(R 2 mT)2, with a 3> 1. In this case, 
zV'(z) =M--Lr 'dy ^jfl!f = A/" + aLia . 



(4.44) 



Based on what we have learnt hitherto, we expect that in this case there is a 
confined phase for small enough fi and then a transition to a deconfined phase at 
fi ~ m. In this case for larger fi the deconfined phase persists since the cut in ( |4.44j ) 
does not allow the end-points of the open distribution to close up again. So the 
validity of the continuum approximation implies a T which is large enough so that 
there is no series of transitions. The confined phase is simple to analyse using the 
techniques in earlier sections. One finds that M = and the density is 



CT 



= -- - Li |(-^) • 
The contour C is determined by inverting 



(4.45) 



is -<xLi5(-£z) 

e = ze 2 



(4.46) 



The phase persists up to a value of fi such that the zero of g(z) at, say z = z 
(Imz — and z < 0), lies on the contour C, which means that z and £ are roots of 



1 — aLis (— £z ) 



n j -<xLi 5 (-£zo) 1 

and z e ? = — 1 . 



(4.47) 



Since a ^> 1, the solution is approximately z$ = —e and £ = (ecr) 1 and so the 
transition occurs at 



// 111 - T log( err) 111 - T 1 + - log(-R 2 mT) + - log 2tc + log 



(4.48) 



The deconfined phase in this case is described by a resolvent 



^ Jo 1 + £e-^ 



V 7 ^ - I)(z - I*) 



(4.49) 



In principle, one can find z numerically, however, we don't pursue this analysis here. 
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5. The Large N Theory in the (//, T) Plane 



When we move out of the regime of low temperature, so that T is no longer <C Br 1 , 
into the whole (p, T) plane, the factors z b {nf3/ B) in ( j2.20| ) cannot be ignored. Hence, 
the measure on P can no longer be approximated by that of a unitary matrix. In this 
case, we can in principle solve the model in the confining phase, where the eigenvalues 
lie on a closed contour by introducing the Laurent expansion 

oo 

p(z) = P^' 1 ' 

n=— oo 

where po = 1 follows from the normalization condition. In terms of the components, 
the action of the Polyakov line with Nf fundamental fermions is 

oo 1 

S(p n ) = N 2 ~ [(1 - z b {nf3/B)) Pn p_ n 

n=i n (5.2) 

+ ^(-l) n z f (n/3/B,mB)(e n ^p n + e""^) • 
The action is quadratic in the p n and so the saddle-point is easily found to be at 



z f(\n 


/3/B,mB) 


1 - z b (\n 





Pn = ^{-\)™e^ r :L* al J > ( 5 - 3 ) 



for n 7^ and po = 1. Note that the SU(N) condition ( [4.10 ) is automatically 



satisfied. The boundary of the confining phase is then obtained by solving for the 
condition that p(z ) = for zq lying on the contour C which itself is determined from 



When m ^> B 1 , and in the intermediate regime 

1 i 

— > T > (5.4) 

R mB 2 V ; 



we have already determined that there is a Gross- Witten transition ( 4.48 ). At p = 
we know that, for m > |, the fermionic modes are decoupled. So we have effec- 
tively the pure gauge theory which is known to have a confinement/deconfinement 
transition driven by the modes p±i at a temperature where z^fi/B) = 1, giving 
T ~ 0.759i? _1 [3]. This transition is known to be first order even when higher order 
corrections in the coupling are considered; in fact the first correction to the Gaussian 
approximation requires a 3-loop computation [39]. For finite fermion mass the tran- 
sition occurs below 0.759-R -1 . If we make the approximation of ignoring the higher 
modes p n , \n\ > 1, then the transition happens when p\ = p_i = 1/2 such that 

1 Nf 1 

-z b (P/B) + jfzfW/B, mB) = -, (5.5) 
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Figure 15: A schematic picture of the phase structure of theory with m S> 
structure of the energy levels is only apparent when T <C which we have indicated by 
the sawtooth pattern. 



and is a third order Gross- Witten transition. The case with m = was considered 
in [40]. 

The Gross- Witten transition at \i = extends out into the (/i, T) plane and joins 
smoothly with the transition line that comes from the first energy level. Although the 
line of transitions can be determined numerically a good approximation is obtained 
by keeping only the first 2 modes p±\. A schematic picture of the phase diagram 



for the large mass regime is shown in Figure [15|. We have also indicated on this 
diagram the discrete sawtooth structure which can be seen at very low temperatures 
T < In this regime the peaks rise until /i approaches ~ 2m and then they 

decrease slowly. 

When m is reduced, the transition at fi = occurs at a lower temperature and 



a schematic picture of the phase diagram appears in Figure [16|. Now the sawtooth 
structure extends all the way to small \x. 



6. Comparison of N = 3 and oo 



It is meaningful to compare the results we have obtained using the two different 
calculational techniques presented in this paper. This is useful not only towards 
checking our results, but also towards learning about what differences might arise 
when considering small vs. large N, and under which conditions there is agreement. 
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Figure 16: A schematic picture of the phase structure of theory with m = based on 
numerical approximations. The deconfined phases touch the T = axis at the points 
H = (i + ^R- 1 , £= 1,2.... 
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Figure 17: Comparison of N = 3 and N = oo QCD at the first level transition (/ = 1) 



with mR = 0, o\ 



(Left): Fermion number. (Right): Polyakov lines &\ and 



We consider a single level transition, i = 1, and <r = |, which is the value used in 
the N = 3 calculations. In Figure [T7| (Left) we show an overlay of the N = 3 and oo 
results for the fermion number normalized by iV 2 . The transitions between the small 
£ confined phase and the deconfined phase, and also between the deconfined phase 
and the large £ confined phase, are clearly signified by discontinuities for N = oo, 
but they are smoothed out for N = 3. Taking lower temperatures does not sharpen 
the N = 3 result since we are plotting as a function of log£ rather than [i. It may be 
that the transitions are smoother for iV = 3 as a result of working in a small spatial 
volume. The important feature of this comparison is that the iV = 3 calculation 
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Figure 18: The points in each figure represent the distribution values Zi = e l( ^ >i of the 
angle expectation values fa of the Polyakov line for N = 3. Four values of £ in each of the 
regions £ < £1, £1 < £ < £2, and £ > £2 are considered. The contours connecting the points 
are interpolations between them and show how the distribution moves with £ (in each plot 
as £ increases the contours change from red— > yellow — )>green— )-blue) . This is in qualitative 
agreement with the large N distributions. 



suggests possible transitions and the N = 00 calculation confirms them. 

In Figure [17| (Right) we show the Polyakov lines and for N = 3 and 00. 
Again the discontinuities which occur at the transitions in the N = 00 case are lost 
for iV = 3. In addition, since the difference \& x - @>*_^ N=co) > \&> 1 - ^'* 1 | (iV=3) , 
and given our observation that the severity of the sign problem increases with 
\&*\ — then it appears that the sign problem would always be more severe 

in the large N limit. 

It is also possible to solve for the distribution zi = e 1 ^ from the expectation 
values fa of the angles of the Polyakov line for N = 3, for the purpose of comparing 
with N = 00. To do this we solve the characteristic equation 

( e *i _ z )( e <fc _ z )( e -^ +9 ^ -z) = 0. (6.1) 

This result can be expanded and rewritten in terms of the expectation values = 
e i6i _|_ e ie 2 _j_ e -i(6i+e 2 )^ anc j op_ x = e -ifli _|_ e -ie 2 _|_ e m+e 2 )^ determined earlier, to give 

z 3 _ z 2^> 1 + z 0>_ 1 -i = o. (6.2) 

We find the three solutions for z (using our results for £P\ and at the I = 1 level 
transition) for four different values of £ in each of the regions £ < £1, £1 < £ < £2, 
and £ > £2- The results are given by the points in Figure |18[ The contours are 
interpolations of the points and give a rough idea of how the distribution changes 



with £. Comparing with the large N results for the distribution, given in Figure 12 



we see that the direction of movement of the distribution matches up in all three 
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regions: In the complex z-plane, for £ < £1 the contour passes from the unit circle to 
the left with increasing £, for £1 < £ < £ 2 it passes back to the right passing by the 
unit circle, and for £ > £ 2 the distribution moves left again back to the unit circle. 

It is also worth noting that there is agreement between the large N finite temper- 
ature phase diagrams in Figures |T3| and [Tj], and the low temperature results for the 



Polyakov line for N = 3 in Figures |2|, and |8| (Left). In the case of a massless quark, 
and at low temperatures, we see a widening of the deconfined phases for increasing 
jj,. For a massive quark, at low but non-zero temperature, there is a brief respite 
from the oscillating phases after the onset transition to the first deconfined phase at 
fi = m, but taking \i large enough causes the oscillations to return. 



7. Conclusions and Outlook 



We have performed a one-loop analytical derivation of the phase diagram of QCD 
as a function of temperature and chemical potential, in the large N and Nj limit, 
on S 1 x S 3 , and supported the low temperature results with numerical calculations 
for N = 3. In the case of massless quarks, and considering the low temperature 
limit, we observe a series of confinement-deconfinement transitions as a function of 
the chemical potential. In the large N limit the phases are characterized by the 
distribution of the Polyakov line eigenvalues in the complex plane which can be 
obtained using matrix model techniques that have been generalized for a complex 
action. In the large quark mass limit we observe the "Silver Blaze" feature in that 
bulk observables are roughly zero until the onset transition to the deconfined phase 
which occurs at /i = m. From here there is a brief continuum-like behavior in 
that the observables appear smooth after the transition. This behavior continues 
until fi is sufficiently larger than m that the levels spread out and the confinement- 
deconfinement phase oscillations return. 

From a physical standpoint, a remarkable aspect of our results is the correlation 
between deconfinement and the existence of partially filled quark shells in a box. 
From the behaviour of the fermion number at low temperatures we observe that 
each confinement-deconfinement transition in the Polyakov line is associated with 
a level transition. All energy levels below the Fermi energy, given by the chemical 
potential, are filled, and levels above the Fermi energy are empty. There are clear 
analogies with partially-filled bands in condensed matter physics; here we would say 
that the Fermi energy in such a system falls within a band, resulting in a conducting 
ground state. A more accurate analogy, since it involves a finite number of particles, 
would be partially-filled shells in nuclear physics. The lesson we draw is that de- 



- 39 - 



confinement appears to require a non-zero density of gapless states. We also learn 
that these states are either particle-like or hole-like depending on which edge of the 
band is closer. Because particles and holes carry conjugate representations of the 
gauge group the resulting physics is distinct, as revealed by the differing behaviours 
of &\ and The non-monotonic behaviour of ^(/i) has been observed in lattice 

simulations of QCD with gauge group SU(2) near its saturation density (i.e. 2N c Nf 
quarks per lattice site) [20]. 

There are a number of interesting generalizations that might be made. One 
could consider more quark flavours with different masses rrif, perhaps coupled to 
different chemical potentials \x$. A detailed multidimensional phase diagram could 
be calculated as a function of the quark masses and their chemical potentials. There 
may be a sign of color superconducting phases from the configurations of the gauge 
field. It might be interesting to consider different manifolds as well. In addition, a 
numerical computation of the line of transitions in the (/i, T) plane can be performed 
from the results for the large N calculation to connect the transition points on the 
T and \x axes for large quark mass, and the curvature of this line can be determined. 
In the event that lattice simulations of QCD become possible at moderate chemical 
potential in the low temperature limit, it would be interesting to see if the low 
temperature confinement-deconfinement phase oscillations we have observed as a 
function of /i are present at strong coupling and/or large volume. We are currently 
investigating this possibility in simulations of QCD with N = 2. 

Our analysis has been limited to the Gaussian approximation and it is clearly 
important to consider systematically the effect of higher orders in the gauge coupling. 
This kind of higher order analysis was qualitatively undertaken in the thermal M = 4 
theory in [3]. It turns out that in this case the fate of the confinement /deconfinement 
transition depends critically on the sign of the coefficient of a 3-loop term in the ex- 
pansion of the effective action. Whilst this 3-loop calculation has not been performed 
for the Af = 4 gauge theory it has for the pure gauge theory in [39]. This tour de 
force calculation proves that the transition survives as a first order transition. The 
effect of higher orders in the coupling for the theory with fundamental flavours (with 
H — 0) have been considered qualitatively in [41,42]. This latter work suggests in 
the theory with /i = that the confinement /deconfinement transition disappears 
and becomes a cross-over as is increased beyond a certain critical value. It will 
be interesting to see how this is altered in the presence of a chemical potential. It 
will be interesting to see whether the phase structure that we find can be related to 
that of a theory with an AdS/CFT-type gravity dual. For instance, will the infinite 
sequence of Gross- Witten transitions that we see be seen in the dual gravitational 
description? 
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